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and 

Laboratory  for  Surface  Science  and  Technology 

University  of  Maine 
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Abstract 

'We  use  a  modified  Kadanoff's  variational  method  to  calculate  the  phase 
diagram  of  an  Ising  model  with  random  vacancies  on  one  of  two  interpenetrating 
sublattices  of  the  isotropic  square  (SQ)  and  body-centered  cubic  (BCC)  lattices. 

We  find  second  order  phase  transitions  only  for  T>0.  The  transition  temperature 
to  very  good  approximation  decreases  linearly  with  impurity  (i.e.  vacancy)  concen¬ 
tration  at  small  concentration.  This  agrees  with  the  linear  decrease  observed  in 
other  systems.  A  plausible  explanation  of  the  absence  of  first  order  transitions 
for  T>0  is  given. 

The  relation  of  our  models  to  certain  percolation  problems  is  similar  to  that 
of  some  spin  models  studied  by  Syozi .  Thus  our  results  allow  an  estimate  of 
critical  probabilities  for  percolation.  . _  / 

i 

1 .  Introduction 

In  recent  years,  much  attention  has  been  paid  to  extensions  of  the  ordinary 

spin  1/2  Ising  model,  such  as  the  Potts  model  and  spin  one  models  An 
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example  of  the  latter  is  the  Blume-Capel  model  ’  represented  by  the  Hamiltonian: 

H  =  -J 7.  S.  S.  -  D:  (1-S.2)  -  B?:S .  (1) 

1  J  i  1  i  1 

<nn> 

where  s^,  Sj  =  1,  0,  -1,  £nn>  is  a  sum  over  nearest  neighbor  s.,  s^  on  the  lattice, 
and  D  is  the  zero-field  splitting  parameter.  Eq.  (1)  may  be  considered  as  an  Ising 


model  with  annealed  impurities  (or  random  vacancies)  represented  by  the  zero 
component  of  spins  s^.  The  thermal  average  impurity  concentration  is  controlled 
by  the  parameter  D.  As  D/0  the  system  reduces  to  the  pure  Ising  model.  Mean 
field  Monte-Carlo®,  and  renormalization  group  (R  G)  studies  of  Eq.  (1) 
indicate  that  the  order-disorder  phase  transitions  of  this  model  consist  of  a  line 
which  is  first  order  at  low  temperatures  and  second  order  at  high  temperatures.  It 
is  generally  agreed  that  two  regions  meet  at  a  tricritical  point;  however,  the 
location  of  this  point  is  not  exactly  known. 

1 2-1 5 

In  this  paper,  we  apply  a  modified  Kadanoff's  Variational  method  (MKVM) 
to  calculate  the  phase  diagram  of  a  related  model  given  by: 

H  =  -jEa.  s.  -  Dr  (1-s.2)  -Be  o.  -  BE  s.,  (2) 

<nn>  3  j  3  i  1  j  3 

where  a.  =  +  1,  s.  =  1,  0,  -1,  <nn>  denotes  nearest  neighbors  as  in  Eq.  (11),  J>0 
and  D  is  a  zero-field  splitting  parameter  for  s..  We  consider  Eq.  (2)  for  $Q  and 

J 

BCC  lattices  in  two  and  three  dimensions  respectively.  We  give  results  for  the 

B=0  case  only.  The  spins  s.  are  located  on  one  of  two  interpenetrating  sublattices, 

J 

the  spins  on  the  other.  The  two  dimensional  case  is  illustrated  in  Fig.  1.  In 
three  dimensions  the  spins  oi  sit  on  a  simple  cubic  lattice  with  a  spin  s..  at  the 
center  of  each  unit  cell.  Eq  (2)  may  be  considered  as  an  Ising  model  with  annealed 
impurities  or  vacancies  present  on  any  site  j  of  one  sublattice  when  s-  =  0.  The 

J 

thermal  average  of  impurities  is  controlled  by  D  as  in  Eq.  (  *  D/J->--«>,  Eq  (2) 

also  reduces  to  a  pure  Ising  model.  In  sect.  II,  we  calculate  phase  diagrams  of 
our  model  and  find  that  the  order-disorder  transitions  of  Eq  (2)  for  T>0  consist  of 
a  line  of  second  order  transitions  only.  The  first  order  line  found  in  the  phase 
boundary  calculated  from  Eq.  (1)  is  no  longer  present.  A  plausible  argument  for 
this  difference  is  given  .  In  sect.  Ill,  we  consider  our  models  at  T->0.  Our  models 
are  closely  related  to  some  systems  studied  by  Syozi^®.  Each  of  these  spin  models 
has  a  corresponding  percolation  problem.  The  spin  system  phase  boundary  for  the 
order-disorder  transition  at  T->0  gives  an  accurate  estimate  of  the  critical  proba¬ 
bility  for  the  corresponding  percolation  problem.  Thus,  the  critical  percolation 
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probabi 1 i ti es  for  our  models  may  be  estimated  from  our  calculated  phase  boundaries 
at  T->0. 


II.  Method  of  Calculation  and  Results 


The  method  of  calculation  is  essentially  the  same  as  that  reported  in  previous 


13-14 


papers  ”  .  We  briefly  review  it  here.  We  consider  Eq  (2)  for  two  dimensional 
square  (SQ)  and  three  dimensional  body-centered  cubic  (BCC)  lattices.  The  Hamil¬ 
tonian  of  Eq  (2)  may  be  rewritten  as: 


H  =  -ZV  , 
R  K 


VR  =  JSR  (o.|  +...+az)  +D  (1-SR  )  +  B  (o-j  +...+o^)  +  BSR,  (4) 

where  R  is  a  unit  hypercube  of  the  a  spin  lattice  with  £  (=2d,  d  being  the  space 
dimension)  a  spins  at  the  corners  and  one  s  spin  in  the  center.  This  is  illustrated 
in  Fig.  1  for  the  two  dimensional  case.  Performing  an  exact  decimation  calculation 
for  the  partition  function  of  Eq.  (3),  which  sums  over  the  central  s  spin  in  each 
cell,  we  obtain  an  effective  Hamiltonian  for  the  remaining  o  spins.  The  resulting 
unit  cell  potential  is: 

Veff  =  eB^al  +-**+az)/z-  +  3D  (5) 

+  In  {1  +  2  exp  (-pD)  cosh  (SB  +  00  (a^  +...+az)}, 
where  0  =  In  two  dimensions,  for  the  case  B  =  0,  which  we  examine  below,  Eq  (5) 
is  equivalent  to  a  o  spin  Hamiltonian  with  nearest  neighbor,  next  nearest  neighbor, 
and  four  spins  (around  a  unit  hypercube)  interactions.  In  three  dimensions,  there 
are  also  six  and  eight  spin  interactions.  For  given  J,  D  and  0,  we  then  use  Eq  (5) 
as  ^0,  the  initial  cell  potential,  and  carry  out  the  step  by  step  RG  transformation. 

We  calculate  the  free  energy  (f),  internal  energy  (U)  and  spontaneous  magnetization  (M) 
per  spin  (here  M  «<c.  >  +  <Sj>).  The  average  occupation  number  of  the  zero  component 
of  s  spin  <n  >,  ie  the  concentration  of  impurities  or  vacancies,  is  given  by: 


<n0>  =  "  a(eo)  =  \*)  '  si2> 

f  0  I  i  c 
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I 

<n0>  is  calculated  using  the  chain  rule,  by  the  same  method  used  to  calculate  U 
and  M  (Eq  (17)  of  paper  14  with  q  representing  &D).  AS  D/J  +  -*  or  ®,  <nQ>  +  0  or  0.5 
respectively.  From  the  behavior  of  the  coupling  constants  (tending  to  larger  or 

i 

smaller  values),  we  know  whether  the  system  is  in. an  ordered  or  disordered  state. 

From  a  simple  ground  state  energy  argument,  it  is  obvious  that  for  B  =  0  the  order- 
disorder  transitions  occur  only  in  the  region  ^  <  1.  The  MKVM  is  known  to 

give  accurate  values  for  the  transition  temperature  over  a  wide  range  of  parameter 
values. 

The  calculated  phase  diagrams  (for  B  =  0)  are  shown  in  Fig.  2,  where  Tq  is  the 
transition  temperature  of  the  pure  Ising  model  calculated  with  the  MKVM  method, 
which  gives  =  0.4576223  and  0.1621341  for  the  SQ  and  BCC  lattices  respectively. 

0  "I  -7  -Ig  J 

For  comparison,  the  exact  or  series  expansion  values  of  tw-  for  SQ  and  BCC 

K  o 

lattices  are  0.44068679...  and  0.15740  +  0.0003  respectively.  The  part  of  the 
diagram  under  the  curve  is  in  an  ordered  state  and  above  the  curve  in  a  disordered 
state.  We  calculated  U,  M  and  <nQ>  near  the  phase  boundary  and  found  that  there  is 
no  discontinuity  across  the  phase  boundary  even  as  1  and  +  0*  F°r  example, 
for  SQ  lattices  at  =  0.9925,  we  found  that  All  =  1 0 ~ ^ &0 5  A<nQ>  =  10"^  and 

AM  -  10  where  AQ  =  Q  (T+fe)  -  Q  (Tc-t)  with  Q  being  U,  <nQ>  or  M,  Tc  being  the 

D 

phase  transition  temperature  for  given  and  L  being kvery  small  number.  The  aQ 
values  just  mentioned  are  of  the  same  order  of  magnitude  as  numerical  errors 
arising  from  truncation  and  the  approximate  location  of  Tc,  etc.  We  may  thus  con¬ 
clude  that  the  transition  is  always  second  order  rather  than  first  order  for  T>0. 

In  Fig.  3,  we  have  converted  the  phase  diagram  of  Fig.  2  to  Tc/T0  versus  <no>.  It 
is  interesting  to  note  that  as  <nQ>  approaches  the  critical  concentration  <nQ>c, 
above  which  there  is  no  phase  transition  at  all,  Tc  decreases  smoothly  to  zero  at 
<no>c  *  0*26  in  the  SQ  lattice  while  Tc  jumps  rapidly  to  zero  near  <nQ>c  =  0.42  in 

the  BCC  lattice.  In  the  later  case  U  and  M  still  change  continuously  when  crossing 

6  8 

the  phase  boundary.  Mean  field  and  Monte-Carlo  studies  of  the  phase  diagram  for 
the  Blume-Capel  Model  (on  a  SQ  lattice)  indicate  that  at  T  =  0i<nQ>c  for  this  model 
is  1.0. 
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The  part  of  Fig.  3  for  small  <nQ>  is  enlarged  and  reproduced  in  Fig.  4,  which 
indicates  that,  to  very  good  approximation,  Tc  decreases  linearly  with  <nQ>  for 

small  <nQ>.  The  slopes  of  the  straight  lines  in  Fig.  4  are  1.683  +  0.002  and 

/ 

1.116  +  0.001  for  SQ  and  BCC  lattices,  respectively.  The  phase  transition  temper¬ 
atures  of  the  Ising  model  with  quenched  ^  or  annealed  ^  bond  impurities  have  also 
been  studied  by  series  expansion  method  and  are  found  to  decrease  linearly  with 
small  impurity  concentration.  It  has  also  been  observed  experimentally  that  in 
nickel-rich,  nickel-titanium  alloys  the  ferromagnetic  Curie  temperature  decreases, 
to  very  good  approximation,  linearly  with  titanium  concentration  x,  for  x  -  1.5%. 

The  approximately  linear  decreast  of  T  with  small  impurity  concentration  is  perhaps 

c 

a  rather  general  phenomena. 

Now,  we  give  a  plausible  explanation  for  the  absence  of  a  first  order  phase 
transition  in  the  model  of  Eq.  (2)  for  T>0.  A  system  at  a  second  order  phase 
transition  is  characterized  by  an  infinite  correlation  length  £,  ie  a  perturbation 
at  one  point  can  propagate  to  points  an  infinite  distance  away.  In  the  model  of 
Eq.  (2),  the  impurities  (zero  component  of  s.)  appear  only  on  one  sublattice.  In 

J 

the  region  of  interest,  <  1 ,  the  spins  o  on  the  other  sublattice  always  have  non¬ 
zero  effective  coupling  (see  Eq.  (5)).  Thus  the  correlation  length  can  easily  extend 
to  infinity  at  Tc,  in  contrast  to  the  model  of  Eq.  (1),  where  impurities  can  appear 
in  both  sublattices  and  thus  more  easily  destroy  the  propagation  of  a  perturbation 
from  one  point  to  another  far  away  point. 

Vie  can  also  understand  the  absence  of  first  order  transitions  as  follows.  For 
a  square  lattice  one  expects  a  first  order  transition  for  <  1,  at  low  T  values, 
if  it  exists  anywhere.  It  follows  from  Eq.  (5)  that  the  (exact)  decimation  of  the 
spins  $.j  generates  (for  B  =  0)  effective  nearest  neighbor,  next  nearest  neighbor, 
and  four  spin  interactions  that  we  denote  J  ,  J  ,  and  J.  respectively.  By 
permutation  symmetry  one  has  J  „  =  J  „  .  In  general,  J  >0  for  ferromagnetic 
J  (J  >  0).  In  the  region  of  possible  first  order  transitions  one  can  show  that 
J4  >  0  as  well.  Thus  when  the  vacancies  are  summed  away  all  interactions  promote 
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ferromagnetic  ordering  of  the  a  spins  in  this  region,  and  a  second  order  transition 
(of  simple  Ising  type)  is  to  be  expected,  as  we  find  in  our  MKVM  calculation.  The 
conclusion  of  the  absence  of  first  order  transitions  on  a  BCC  lattice  for  gj  <  1 

i 

can  also  be  reached  by  a  similar  argument. 

III.  Low  Temperature  Limits 

Syozi  has  calculated  the  exact  phase  diagrams  of  Ising  models  on  decorated 

flew«L  1 

square  (SQ),  honeycomb  ( HC  )  and^triangular  (T’T)  lattices  with  spins  s  =  ^  on 

each  lattice  site  and  spins  s  =  1  on  each  nearest  neighbor  bond.  He  also  calculated 

A 

phase  diagrams  for  Ising  models  onAsemi-dil(  te  honeycomb  lattice  where  the  spins 
s  *  1  and  s  =  1  spins  are  on  two  equivalent  interpenetrating  sublattices  as  in  the 
present  case.  The  phase  boundaries  of  Syozi's  models  at  T-+0  give  the  critical 
concentration  Pc  of  spins  (1-P  is  the  vacancy  concentraiton  )  on  sites  of  s  =  1 
spins  above  which  ferromagnetism  can  occur.  It  turns  out  that  Pc  calculated  in  this 
way  is  exactly  equal  or  very  close  to  the  critical  probability  of  the  corresponding 

pp 

bond  (for  decorated  lattices)  or  sublattice-site  (for.  semi-dilute  HC  lattice  ) 
percolation  problems.  For  the  Syozi  model: 

Pc  (decorated  SQ)  =  0.5 

P  (decorated  pTl  =  1-  P  (decorated  HC)  =  ^ =  0.3522  ... 

C  C  lo 

Pc  (semi-dilute  He)  =  0.5. 

For  the  percolation  problem  22,  23. 

Pc  (bond,  SQ)  =  0.5 

Pc  (bond,  PT)  *  1-  Pc  (bond,  HC)  =  2  sin  =  0.349... 

Pc  (sublattice  site,  HC)  =  0.5  (8) 

It  is  obvious  that  our  model  of  Eq  (2)  on  the  SQ  and  BCC  lattices  is  just  an 
Ising  model  on  a  seini-dilute  SQ  or  BCC  lattice  and  its  low  temperature  behavior  is 
similar  to  the  Syozi  models.  Thus,  the  critical  probabilities  PG  of  the  sublattice- 
site  percolation  problem  on  the  SQ  and  BCC  lattices  can  be  estimated  from  <no>c  of 


Applying  Eq.  (9)  to  <n0>c  of  Fig.  3,  we  have  Pc  55  0.48,  0.16  for  the  SQ  and  BCC 
lattices  respectively. 
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Fig.  1.  Lattice  of  Eq.  (2)  in  two  dimensions.  Dots  (.)  spins  o.. ,  crosses  (x) 

spins  s..  The  dotted  lines  enclose  a  unit  hypercube  of  the  o  spin  lattice. 
J  , 

Fig.  2.  Transition  temperature  T  vs.  D/Jz  obtained  for  a  two  dimensional  SQ  (+) 

c 

lattice  and  a  three  dimensional  BCC  (a)  lattice. 

Fig.  3.  Transition  temperature  vs.  impurity  concentration  <nQ>  for  a  SQ  (+)  and  a 
BCC  (a)  lattice. 

Fig.  4.  Part  of  the  curves  in  Fig.  3  for  <nQ>  <2.6%.  The  solid  curves  are  obtained 
by  a  linear  least  square  fit  to  the  data  points. 
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